Single-qubit unitary gates by graph scattering 
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We consider the effects of plane-wave states scattering off finite graphs, as an approach to im- 
plementing single-qubit unitary operations within the continuous-time quantum walk framework of 
universal quantum computation. Four semi-infinite tails are attached at arbitrary points of a given 
graph, representing the input and output registers of a single qubit. For a range of momentum eigen- 
states, we enumerate all of the graphs with up to n = 9 vertices for which the scattering implements 
a single-qubit gate. As n increases, the number of new unitary operations increases exponentially, 
and for n > 6 the majority correspond to rotations about axes distributed roughly uniformly across 
the Bloch sphere. Rotations by both rational and irrational multiples of n are found. 



I. INTRODUCTION 

Quantum walks are the quantum mechanical analogs 
of random walks in classical systems, and like their clas- 
sical counterparts they have provided a useful framework 
for the construction of efficient quantum algorithms [l]- 
In addition to the development of quantum walk 
approaches to known quantum algorithms such as the 
Grover search Q, triangle finding @, and element dis- 
tinctness 0, several quantum algorithms were obtained 
using quantum walks that were previously not found 
within the standard quantum circuit or adversary mod- 
els, such as the traversal of glued-tree graphs with ran- 
domness dl , and the solution of NAND trees [1, , min- 



max trees and general boolean functions [ij 



More recently it has been shown that quantum walks 
are able to perform arbitrary quantum computational 
tasks ^M^- 

In the continuous-time approach to universal quantum 
walks computational basis states are represented 

by long linear graphs (tails) which support plane wave 
modes characterized by momentum eigenstates. A single 
quantum walker is initialized in a momentum eigenstate 
of a single semi-infinite line out of 2^ such tails repre- 
senting N qubits. Quantum operations on the encoded 
information are carried out by interspersing the tails with 
small graphs. The scattering of the plane waves off these 
graphs mixes the amplitudes on different tails and adds 
possible phase shifts. 

In Ref. a small number of such graphs were iden- 
tified that could generate a universal set of A^-qubit 
unitary gates for a fixed momentum k ~ — 7r/4 (all 
length units are suppressed for convenience). The two 
graphs needed for single-qubit operations are shown in 
Fig. [T] The first of these, depicted in Fig. [Tal is a 
disconnected seven- vertex graph that yields the unitary 
Rz{~n/4) ~ diag(l, e~'^/*) = T neglecting unimpor- 
tant overall phases {Ra{6) denotes a rotation about axis 
a by an angle 9, and X, Y , and Z are the Pauli matri- 
ces). The upper subgraph is a two-site graph, so with 



the tails encoding the |l)input and 



l)output attached it 



has the effect of changing the length of the full tail by 
one lattice spacing. For plane waves of the form e''^^, 
this corresponds to a scattering phase shift of — 7r/4 on 
the |1) state. The lower subgraph, attached to the |0) 
input and output tails on the same vertex, produces no 
scattering phase shift and one thus obtains the desired 
rotation. This subgraph furthermore has the desirable 
property of also having an effective length of one lattice 
spacing, so that the registers |0) and |1) maintain their 
spatial coherence after the scattering event. Under this 
scheme the various paths through a graph be of equal 
length, and there must be an identity gate with the same 
effective length for a given momentum to act on all regis- 
ters where no gate is desired. This identity requirement 
is trivial to satisfy in the case of graphs with integral ef- 
fective lengths, but as we show below graphs can have 
lengths that are non-integral, or even irrational. 

The second graph, shown in Fig. Ilb| is a connected 
graph that yields Rx(tt/2) = V, again neglecting unim- 
portant overall phases. The Hadamard gate can be con- 
structed via H = T^VT^, and combinations of H and 
T can produce any single-qubit unitary at fixed accu- 
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racy }1YIJ . While these two graphs are sufficient to gener- 
ate a universal set of single-qubit gates, and a universal 
set of gates when an entangling gate is included , it is 
not clear if there are other graphs that could be chosen 
instead. 

An immediate question is: are there other graph com- 
binations that work equally well for k = — 7r/4, or for 
other values of the momentum? Because the number of 
non-isomorphic graphs grows exponentially in the num- 
ber of vertices n, one might expect a plethora of choices 
which would be useful in the design of quantum algo- 
rithms based on quantum walks. Yet it is also conceiv- 
able that larger graphs will simply reproduce the uni- 
taries found for smaller graphs, or that as n increases 
the requirement that there be perfect transmission from 
each of the input tails to the combination of the two out- 
puts, with equal effective lengths on all four paths, will 
become increasingly difficult to satisfy. 

In the present work, we explore these issues numer- 
ically by calculating the single-qubit unitary operators 
resulting from two input tails scattering to two output 
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FIG. 1. The graphs identified in Ref. that generate a uni- 
versal set of one-qubit unitary gates, with vertices denoted 
by circles and edges by solid lines. Open circles indicate 
that both an input and an output tail attach at that ver- 
tex, whereas circles open only to the left or right indicate at- 
tach me nt of only input or output tails, respectively. Graphs 
(a) and (b) generate the rotations J?z(7r/4) and 7?x(7i"/2), re- 
spectively, for k = — 7r/4. The tails, depicted by dashed edges 
and ellipses to denote their continuation, are included to in- 
dicate the attachment vertices. These graphs appear in the 
supplemental material as IDs 238 and 309. 



tails, off every possible graph with up to nine vertices. 
Momentum eigenstates k — {p/q)Ti where q G {2,3,4,5} 
and Q < p < q for integer p are considered; we study 
only positive values of k for two reasons. Mathemati- 
cally there is a symmetry between positive and negative 
momenta such that negating k conjugates the reflection 
and transmission coefhcients; physically if one considers 
propagating wavepackets with tightly peaked momenta 
instead of un-normalizable plane waves, the negative mo- 
menta result in the initial state traveling away from the 
scattering graph instead of toward it, under the defini- 
tions we use. Note that in Ref. Q the graphs of Fig. [T] 
are utilized at momentum k — — 7r/4 because the Hamil- 
tonian employed therein is equal to the adjacency matrix 
of the graph, rather than its negative. We find that the 
number of graphs meeting the requirements for imple- 
menting single-qubit rotations does indeed grow expo- 
nentially with n. The most prolific momentum we inves- 
tigate IS k = 27r/3, which results in a total of 98 unique 
unitary operations on graphs of nine or fewer vertices. 

There is a supplemental data file associated with this 
manuscript, available on the arXiv pre-print server [l8| . 
and our main focus here is the analysis of those data. 
This manuscript is organized as follows. The mathemat- 
ical formalism underlying the scattering of plane waves 
by graphs is briefly reviewed in Sec. followed by a de- 
scription of the numerical strategy to enumerate graphs. 
The results of the survey are given in Sec. IIIIl and are 
summarized in Sec. IIVI 



II. FORMALISM 

A. Graph Scattering 

The theory underpinning the scattering of plane waves 
by graphs has been discussed in detail by Varbanov and 



Brun [l9[ , so it will only be sketched here. One considers 
a finite graph G = {E, V} with V a set of n vertices and 
E C V X V a. set of edges defined by the graph adja- 
cency matrix A"-^ . Each vertex v £ V has attached to 
it My semi-infinite tails (My may vanish for as many as 
n — 1 of the vertices of G), with the adjacency matrix of 
tail m on vertex v defined as [^Jm] ■■ = j+i + 
where i,j G {1,2,3,...}. Vertices on the tails are de- 
noted \jvm), with the attachment vertex on the graph 
labeled \Oy„i) where the index m is kept only to match 
the tail labeling, i.e. \Oym) = b)- The adjacency matrix 
associated with the connection of the tails to the graph is 
-4.™ = |0„m)(l„™| + \ lvm){Ovm\, and the total adjacency 
matrix for the scattering system as a whole is 



A = A'- 
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To describe single-qubit unitary transformations we are 
interested only in those cases with four tails — two input 
and two output — attached to G. That is, we require 



Y^My^i. 



(2) 



For the moment though, suppose that the graph G 
consists of a single point. In this case, two semi-infinite 
tails corresponding to the input and output states fuse 
into a single infinite line. The (unnormalized) eigenstates 
of the graph Hamiltonian H = — A are plane waves 
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(3) 



with eigenvalue = — 2cos(fc), where k G [— 7r,7r) is 
the momentum quantum number. One can envisage the 
input as a wavepacket strongly peaked at momentum k 
scattering off a single point, resulting in an identical out- 
put wavepacket. 

Now consider an arbitrary number of tails connected to 
an arbitrary finite graph G, again defined by adjacency 
matrix A and corresponding Hamiltonian H ~ —A. One 
again assumes that the input and output states are plane- 
wave eigenstates of their respective semi-infinite lines, 
with eigenvalue 2cos(fc); that is, the scattering has not 
changed the value of fc, equivalent to the conservation 
of momentum and energy. In principle the graph could 
mix momenta on different input and output tails; this 
possibility is not included in the analysis, because it 
would imply that the quantum information encoded in 
the wavepackets would travel down the tail at different 
rates. 

Consider a momentum eigenstate \k,vm) with incom- 
ing component only on tail m of vertex v, and possible 
outgoing components on every attached tail including the 
incoming one. This leads to transmission and reflection 

respectively, defined by the 
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expressions 



{3vm\k,vm) 
{jv'r,i'\k,vm) 



v'm' .vm 



(fc)e''=J'. 



ikj 



(4) 
(5) 



The first and second terms in Eq. (|3]) represent tlie in- 
coming and reflected component of the state respectively. 
The term on the right hand side of Eq. ^ represents the 
component of the incoming state that is transmitted to 
the other tails (for which v' ^ v or m' ^ m ). These are 
Eqs. (4) and (5) of Varbanov and Brun [l9j . 

The momentum eigenstates |fc, vm) can be expressed in 
terms of their components on the graph vertices |fc, vm)'^ 
and the tails: 



\k,vm) = |fc,i;m)^ + ^(e + r„™e''=J')|j„„) 
/ ^ 

~l~ ^ ^ iy' m' .vm ^ ^ ^ '^\Jv'm'}-j (6) 
v' ,ni' j — ^ 

where the prime on the second sum indicates that all 
possible values of v' and m' for which v' ^ v or m' ^ m 
are being summed over. Enforcing the condition that 
H\k,vm) = — 2 cos(fc)|fc, vm) and simplifying, one ob- 
tains that IkjVm)^ must satisfy 

^A^ - 2cos(fc) +e*'=^M„- \v') {v'\^ \k,vmf 

= 2isin(fc) |u) . (7) 

This is Eq. (10) of Varbanov and Brun [l^. The so- 
lution to the scattering problem consists of determining 
the vector \k,vm)'^ and from it obtaining the reflection 
and transmission coefBcients. Because the reflection and 
transmission coefficients at the attachment points must 
be the same as on the tails, Eqs. ([4]) and ([5|) can be ex- 
pressed more conveniently in terms of \k,vm)^ and the 
graph vertex states as 



{v\k, vm) ' = 1 H 
{v'\k, vm)'^ ~ tyi' 



(8) 
(9) 



In order to interpret the scattering event as implement- 
ing a singie-qubit operation, one assumes that the finite 
graph G has four semi-infinite tails attached. Two of 
the tails are considered as inputs and two as outputs. A 
continuous-time quantum walker in a wavepacket state 
with momentum tightly peaked about k propagates to- 
ward G along the two input tails. It scatters from an ini- 
tial superposition supported on these tails to a superpo- 
sition with support on the two output tails. Any relative 
phase shift or change in probability amplitude between 
the two tails corresponds to a unitary transformation of 
the input computational state. In general, a momentum 
eigenstate incoming on the first input, vm = Oin, can 
be thought of as mapping incoming to outgoing states 



according to 

\k,Oin) ^ ro.„(fc)|fc,Oi„) +fi,„,o.„(fc)|fc,li 



1 ^out/ ■ 

i=0 



(10a) 



Similarly, a momentum state with incoming component 
on the second input, vm = lin, scatters to outgoing com- 
ponents as 

\k,lin) ^ ri,„(/c)|fc,li„)+io.„,i.„(^)|fc,Oin) 
1 

+ Y.*^a..,hn{k)\k,iout), (10b) 

i=0 

SO in order for these scattering processes to correspond 
to a unitary operation mapping |A:,Oin) and Ijn) to 
superpositions of |A:,Oout) and |fc, lout) wc must identify 
combinations of graphs G; tail attachment points |Oin), 
I lin), |Oout), and I lout); and momenta k such that there is 
neither reflection along either input tail nor transmission 
from one input tail to the other. That is, we require 



tn 



0. 



(11) 



When the conditions Eq. pT|) are met, the transforma- 
tions Eq. (fTO]) simplify to 



\k,Oin) ^tOa^t,o,„ik)\k,Oout) +ti„„t.o,„(A:)|A;,lout), 

(12a) 

\k, lin) ^ iOo„t.l,„(fc)|fc,Oout) +tl<,„t,l,„(fc)|fc,lout), 

(12b) 

i.e. an arbitrary input vector a|fc, Oin) + /3\k, Ijn) is trans- 
formed by the operator: 



O 



*0out,0i„ ioout,ii„ 

*lout,Oi„ *lout,lin 



(13) 



We must now show that this operator is unitary. To 
simplify notation, let us label the two input tails by 1 and 
2. and the output tails by 3 and 4. A general scattering 
process on a graph with four attached tails yields an S- 
matrix 



^ ( 



i=i 



Eh 



I 



(14) 



where |Tj) corresponds to the state on tail i, which Var- 
banov and Brun show is unitary [l9j . Our restriction 
that we consider only graphs satisfying Ec^. (jll[l does 
not change the unitarity of S in those cases, and cor- 
responds to ri ~ r2 = ti2 = t2i = 0. We then as- 
sume an initial state with support only on the input tails, 
iV'o) = oi\ti) +/3|r2), with jap + |/3p = 1. The restricted 
S-matrix acts on such a state as 

^l^o) = ("^31 + I3t32) Its) + [aUi + PU2) W), (15) 
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(a) Given two isolated vertices ofT of which to scat- 



FIG. 2. 

ter, the identity and SWAP gates, / and X, are trivial to con- 
struct at any momentum (supplemental material IDs 180 and 
179, respectively, at fc = 7r/4). |(b)| With five vertices it be- 
comes possible to implement gates other than / and X. This 
five-vertex graph implements Rz{—tt/2) and has an effective 
length of 1. (Supplemental material ID 1458, k = n/2.) 



but this is precisely the action of O on a\k, Oin) lin) 
under the identification jr^) -O- \k,i) with 1 o Oin, 2 o 
lin, 3 O Oout, and 4 -H- lout- Since S is unitary in general, 
it remains so when restricted to the cases of interest that 
define O. 



k, {M„}, and A"^ (of order 10^ under the chosen con- 
straints), we solve these equations numerically using the 
GNU Scientific Library [2l| and LAPACK ^ software 
Ubraries. If a particular combination results in reflec- 
tion and transmission coefficients satisfying Eq. pi[) then 
there is zero transmission from the input tail to one of 
the three possible output tails. In this case, we solve the 
equations again with the incoming portion of the mo- 
mentum eigenstate instead on this unused tail. If the 
conditions PT|) are again satisfied, with the same two 
tails supporting the output, then the current combina- 
tion implements a unitary transformation U defined by 
the operator (fT3)) . 

The final step is to determine the effective length 
through the graph from each input to each output. If 
paths of different length exist, then a quantum walker 
initially in a superposition of one wavepacket on either 
input tail would acquire a spatial shift between the two 
output wavepackets after traversing the graph. The ef- 
fective length to output j from input i is defined as [l^ 



dk 



argijout,»i„(fc)- 



(16) 



B. Enumeration of Graphs 

The number of non-isomorphic simple graphs on n 
vertices grows super-exponentially with n, and for each 
graph the number of ways to attach four tails goes as n'^. 
Using the geng {generate graphs) software package, part 
of the popular nauty (no awtomorphisms, yes?) suite 
of graph theoretic tools [2^, we first generate the adja- 
cency matrices of all non-isomorphic simple graphs 
on n vertices for n S {1, 2, 3, . . . , 9}. For each graph 
we then enumerate the set of all combinations of attach- 
ment points for four tails, each of which defines a set 
{Ml, M2, . . . , Mn} satisfying Eq. ©. 

One can see by symmetry arguments that the one- 
vertex graph with four tails attached cannot transmit 
from one of the tails to only two of the remaining three, 
so at least two vertices are required to implement a com- 
putational unitary gate. With two or more vertices avail- 
able, it is trivial to implement the identity gate and the 
SWAP or X gate by constructing a graph with two iso- 
lated vertices and attaching an input and an output tail 
to each one, as exemplified in Fig. I2al If the computa- 
tional labels (i.e. |0) and |1)) agree on each vertex then 
an identity is performed; otherwise the X gate is applied. 
The first unitarics that cannot be trivially produced in 
this manner appear when there are five vertices. An ex- 
ample of such a graph can be found in Fig. I2bl 

For each momentum value k under consideration, for 
each adjacency matrix A'~^ and attachment set {My}, 
Eq. ([7]) produces a set of n linear equations to be solved 
for the coefficients of \k,vm)'^. These in turn yield the 
reflection and transmission coefficients from Eqs. ([5]) and 
Due to the large number of such combinations of 



This can be understood when one considers the case of 
a single input attached to one end of a linear graph of 
L segments, with a single output tail at the other end. 
Clearly the transmission from one end to the other has 
unit magnitude at all momenta, and the phase difference 
between the ends results in a transmission coefficient of 



Equation ([TBI) then yields £{k) = L, as ex- 



pected. We calculate this derivative numerically on each 
candidate graph, using a nine-point stencil, and require 
that the four values £ji{k) agree whenever tj^^^^i^^ ^ 0. 

Our goal is to determine any potential computational 
advantage to increasing the number of vertices in G, so we 
record the current system only if U meets all of the stated 
criteria and has not been found at the same momentum 
and effective length on a smaller graph. If it has been 
previously found, we keep track of how many times it 
has appeared. 



III. RESULTS 

The number of non-isomorphic simple graphs on n ver- 
tices, iV„, does not have a simple closed- form solution. 
Nevertheless, the counts are well documented for many 
values of n (see, for example, Ref. and references 
therein) , and the total number of such graphs on nine or 
fewer vertices is 288 266. For each graph on n vertices, 
the number of distinct ways to attach k tails is 



w„ 



k 
k 



1 



(17) 



and we have chosen to investigate scattering at nine dif- 
ferent momenta with four attached tails. This results 
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FIG. 3. (Top) Total number of graphs resulting in unitary 
operators, as well as the number of distinct unitarics found as 
a function of number of vertices. Here operations are consid- 
ered distinct if they arise from graphs with different effective 
lengths or at different momenta. (Bottom) Number of dis- 
tinct unitaries, as presented in the upper plot, separated by 
momentum. 



n=l 



1262 489148 



(18) 



combinations of momenta, attachment points, and 
graphs to be numerically examined. Of these, it turns 
out that 1 960 316 have the properties required to im- 
plement a computational gate. Figure [3] shows that the 
number of unitary and non-identity operations both in- 
crease super-exponentially in the number of vertices, fol- 
lowing the growth of the total number of graphs, and 
that this is true at each momentum investigated. 

There is significant redundancy within these gates 
however, where many different graphs yield the same 
unitary for a given value of the momentum and effec- 
tive length. This reduces the total number of distinct 
unitary operations to 3380. The adjacency matrices for 
these graphs, the relevant momentum and attachment 
points, and the resulting unitary gates arising, are given 
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FIG. 4. Number of available unitaries provided by increasing 
the number of vertices in the scattering graph, for the five k 
values at which an increase was found. Here we do not con- 
sider unitaries to be distinct if they arise at different effective 
lengths but otherwise perform the same transformation. 



explicitly in the supplemental online information accom- 
panying this work [1^. Of the graphs producing these, 
2496 yield gates other than identities; that said, it is im- 
portant to keep in mind that an identity counterpart is 
required for each non-identity gate in order to ensure uni- 
versality so having a sufficiently large number of identity 
gates is essential. 

Only those graphs which result in unitary operations at 
the same momentum value can be combined under the 
current model to effect a quantum computation, since 
the momentum of the quantum walker is fixed. The ef- 
fective lengths of the constituent graphs are unimpor- 
tant, as long as each graph can be paired with an iden- 
tity operation of the same length. There are 262 graphs 
which do not have commensurate identity graphs, leav- 
ing 3118 candidates for inclusion in computational sets. 
There is further redundancy among these graphs, since 
two graphs performing the same operation with differ- 
ent effective lengths at the same momentum are for our 
purposes equivalent if each has a corresponding identity 
operation. Taking this into account reduces the num- 
ber of potentially useful operations to 284. The number 
of resulting unitaries is shown in Fig. [4] as a function of 
the number of vertices for the momentum values investi- 
gated. There are 16 unitaries that can be produced by 
graphs on five vertices, eight of which are new in that 
they cannot be produced by four or fewer vertices. Sim- 
ilarly, there are 24 new unitaries on six vertices, and 30 
more on seven. These 62 unitary operators are produced 
by only 15 graphs, combined with various tail attachment 
and momentum configurations. One of these graphs has 
two variations of tail attachment points that lead to non- 
isomorphic infinite graphs once the tails are included. 
These 16 graphs are drawn explicitly in Fig. [5] 
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FIG. 5. The graphs on five, six, or seven vertices that im- 
plement single-qubit unitaries which are unavailable on fewer 
vertices. The five- vertex graphs (a) and (b) result in eight dis- 
tinct unitaries under different tail attachment configurat ions 
and different momentum values. The six- vertex graphs (c) 
i)] yiel d a total of 24 unitaries, and the sev en-ve rtex graphs 
lead to 30. Note that while (c) and (d) are isomor- 



Note that while (c) and ^ 
phic, they represent two distinct, non-isomorphic configura- 
tions once tails are attached. 



A. Unitary operations by momentum value 

There are four momentum values, p7r/5 for p = 
1,2,3,4, at which the only available operations are 
the identity and Pauli-X operators; these momenta are 
clearly not useful for our purposes, and are therefore not 
plotted. It is clear that for each momentum value that 
provides a growing number of unitaries, the number of 
unique gates continues the exponentially increasing trend 
as a function of the number of vertices. 

At fc = 7r/4 we find graphs that implement Pauli-X, 
-Y, and -Z gates as well as the rotations i?x(±7r/2) and 
Rz{p7T/i) for p € {— 3, — 2, . . . , 3} (the latter including 
the identity at p = 0). These results contain the graphs 
identified by Childs [IJ] and therefore (re)produce a uni- 
versal set of gates. Additionally, we find graphs produc- 
ing rotations by — tt about those axes in the equatorial 
plane of the Bloch sphere making angles with respect to 
the X axis of pn/S for p £ {1,2,3,5,6,7}. Momentum 



FIG. 6. Distributions of the axes of rotation on the Bloch 
sphere available at momentum values k = n/2, tv/S, and 2-k/3. 
That is, there exists at least one non-trivial rotation R{0^^) 
about each axis defined by the points {9, ip) on the upper 
hemisphere of the Bloch sphere shown here. Symbols are as 
in Fig. H 



k = 37r/4 reproduces this same set of 18 gates. 

When the momentum is k ~ tt/2 wc find graphs ca- 
pable of implementing 52 rotations of the Bloch sphere, 
about 28 non-parallel axes. With momenta A: = 7r/3 and 
fc = 27r/3 we find 90 and 98 rotations about 55 and 59 dis- 
tinct axes, respectively. The available axes are visualized 
for these three momentum values in Fig.[6l The reflection 
symmetries about (f> — and (j) ~ ±7r/2 are due to the 
different configurations of attachment points available for 
each graph. If a given graph performs a single-qubit ro- 
tation about the (9, (f) axis of the Bloch sphere when 
the attachment vertices for (|0)in, |l)in, |0)out, |l)out) are 
|2), |3), |4)), then moving the inputs to the out- 
puts and vice versa by re-attaching the tails in the or- 
der (|3), |4), |2)) results in a rotation about (6', — 0). 
Interchanging the inputs with each other and doing 
the same for the outputs, i.e. re-ordering the tails as 
(|2), |4), |3)), leads to a rotation about (6',7r-(/)). This 
re-ordering is equivalent to conjugating the original gate 
byX. 



B. Rotations by irrational multiples of vr; fractional 
and irrational effective lengths. 

Figure [7] showcases three graphs that produce inter- 
esting and perhaps unexpected results. While there is no 
reason to assume a priori that our procedure will not find 
any rotations by irrational multiples of tt, neither is it in- 
tuitive that this should be the case given that the graphs 
are scattering plane waves with momentum values that 
are rational fractions of tt. Nevertheless, we indeed iden- 
tify over 100 rotations through angles that numerically 
appear to be irrational fractions of tt. These individual 
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(b) 



^5 • • 




FIG. 7. Some graphs that exhibit non-intuitive properties. 
At fc = 7r/3, 1 (a) I performs a Z rotation through an angle of 
arctan (5\/3/llj , followed by a Z gate. Also at fc = 7r/3, 
(b) is an identity gate with an effective length of ^ = 1/2. 
Finally at k = n / 4, (c) performs the basis-changing operation 



VX and |(d) | implements an identity, each with the irrational 
effective length of ^ = 5 — 2^2. (Supplemental material IDs 
802, 488, 327, and 330.) 



results can be checked analytically by solving Eqs. ([7|)- 
(|9]) exactly in the cases of interest; those we checked bore 
out the apparent irrationality indicated by the numeri- 
cal results. Figure [Tal depicts one such case, a graph that 
implements the transformation Zi?2[arctan(5-\/3/ll)] up 
to a global phase at momentum k ~ tt/3. Besides the 
novelty of obtaining irrational multiples of tt under the 
circumstances, these graphs are inherently useful because 
any two rotations about non-parallel axes by irrational 
multiples of tt form a universal set for fixed-precision 
single-qubit quantum computation 
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Similarly, it is not obvious whether any graph should 
have a non-integral effective length. Figure I7bl shows 
a graph on only five vertices that is capable of acting 
as an identity gate with an effective length of £ = 1/2; 
at the momenta in question, no graph exists on fewer 
than five vertices that has a non-integral effective length 
and implements a single-qubit unitary. Finally we see in 
Fig. [7c] a graph that has an irrational effective length of 
£ = 5 — 2\/2. This graph implements the basis-changing 
operation a/X, and remarkably the graph in Fig. I7dl acts 
as an identity gate with this same irrational length. Of 
the 3118 graphs capable of implementing single-qubit 
unitaries with commensurate identities, 2352 of them 



have non-integral effective lengths. Of these, 1042 ap- 
pear numerically to have irrational lengths, including 
the 'longest' graph identified, which has effective length 
e = 350 + 156^5 « 698.826 (ID 2174 in supplemental 
material [3]). This extreme effective length due to a 
comparatively small number of vertices (i.e. 9 ^ 700) is 
another observed phenomenon whose presence is not ini- 
tially obvious. Such a length corresponds to the incoming 
wavepacket's having been localized in the region of the 
graph for a significant duration, and is reminiscent of a 
diverging negative scattering length, approaching unitar- 
ity in traditional quantum scattering theory. Almost 20% 
of the unitaries identified have effective lengths i > 10, 
with greater than 1% having £ > 100. In the same vein, 
it is also noteworthy that no negative effective lengths 
were identified. 



IV. CONCLUSIONS 

Inspired by the proof due to Childs that wavcpackets 
scattering off small simple graphs can perform univer- 
sal quantum computation, we have exhaustively studied 
scattering at nine momenta over the set of graphs on 
fewer than 10 vertices. As the number of vertices in 
the scattering center increases, so does the number of 
distinct single-qubit unitary operations that can be per- 
formed. The number of distinct graphs capable of pro- 
ducing these unitaries grows super-exponcntially, provid- 
ing many methods by which to construct a given opera- 
tor. The promise of investigating graphs on larger num- 
ber of vertices is that a desired operation might be able 
to be implemented by scattering off one or a few graphs, 
where its decomposition into the two transformations in 
a single universal set could require many more. 

Rotations of the Bloch sphere can be performed about 
many axes distributed roughly uniformly across its sur- 
face, by both rational and irrational multiples of tt. This 
results in a large number of distinct sets of operators 
that are universal for quantum computation, which can 
be combined together to minimize the total size of the 
graph used in a computation. Natural extensions to this 
work include increasing the number of vertices under in- 
vestigation, the array of momenta studied, or both; and 
expanding the number of tails attached to each scattering 
graph to find two- or three-qubit gates (or more). 



ACKNOWLEDGMENTS 

We are grateful to Andrew Childs for helpful comments 
in the early stages of this work. The authors acknowl- 
edge funding from Alberta Innovates - Technology Fu- 
tures (MSU), and the Natural Sciences and Engineering 
Research Council of Canada. 



8 



[1] A. Ambainis, Int. J. Quant. Inf. 1, 507 (2003). 

[2] J. Kempe, Contemporary Physics 44, 307 (2003). 

[3] V. Kendon, Math. Struct. Comp. Sci. 17, 1169 (2007). 

[4] V. Kendon, e-print: arXiv: 1107.3795 (2011). 

[5] M. Santha, in Theory and Applications of Models of Com- 
putation (Springer BerUn / Heidelberg, 2008), vol. 4978, 
pp. 31-46. 

[6] F. Magniez, M. Santha, and M. Szegedy, Proceedings of 
the sixteenth annual ACM-SIAM symposium on Discrete 
algorithms (SODA 2005), p. 1109 (2005). 

[7] A. Ambainis, SIAM Journal on Computing 37, 210 
(2007). 

[8] A. M. Childs, R. Cleve, E. Deotto, E. Farhi, S. Gutmann, 
and D. A. Spielman, Proceedings of the 35th ACM Sym- 
posium on Theory of Computing (STOC 2003), p. 59 
(2003). 

[9] E. Farhi, J. Goldstone, S. Gutmann, Theory of Comput- 
ing 4, 169 (2008). 

[10] A. M. Childs, R. Cleve, S. P. Jordan, and D. Yeung, 
Theory of Computing 5, 119 (2009). 

[11] R. Cleve, D. Gavinsky, and D. L. Yonge-Mallo, in Theory 
of Quantum Computation, Communication, and Cryp- 
tography (Springer- Verlag Berlin, Heidelberg, 2008). 



[12] B. W. Reichardt and R. Spalek, Proceedings of the 

40th ACM Symposium on Theory of Computing (STOC 

2008), p. 103 (2008). 
[13] A. Ambainis, A. M. Childs, B. W. Reichardt, R. Spalek, 

and S. Zhang, SIAM J. Comput. 39, 2513 (2010). 
[14] A. M. Childs, Phys. Rev. Lett. 102, 180501 (2009). 
[15] N. B. Lovett, S. Cooper, M. Everitt, M. Trevers, and V. 

Kendon, Phys. Rev. A 81, 042330 (2010). 
[16] M. S. Underwood and D. L. Feder, Phys. Rev. A 82, 

042304 (2010). 

[17] P. Kaye, R. Laflamme, and M. Mosca, An Introduction 
to Quantum Computation, pp. 70-71 (Oxford University 
Press, 2007). 

[18] See Supplemental Material included with this manuscript 
at http://www.arxiv.org for a complete list of unitary 
gates found. 

[19] M. Varbanov and T. A. Brun, Phys. Rev. A 80, 52330 
(2009). 

[20] B. D. McKay, [Fttp : //cs . anu. edu. au/-bdm/nauty[ 

[21] M. Galassi et al., GNU Scientific Library Reference Man- 
ual, http://www.gnu.org/software/gsl 

[22] E. Anderson et al., LAPACK Users' Cmde (Society for 
Industrial and Applied Mathematics, 1999). 

[23] N. J. A. Sloanc, On-Line Encyclopedia of Integer Se- 
quences, http: //oeis . org/A000088j 



